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I.  INTRODUCTION 


There  has  been  increasing  interest  in  the  past  ten  years  in 
the  atmospheric  degradation  of  a  laser  wavefront.  Atmospheric 
turbulence  poses  limitations  on  pointing,  tracking,  imaging 
and  communication.  Researchers  have  attacked  the  problem  from 
two  viewpoints  by  studying  the  turbulence  itself  and  by  making 
direct  measurements  of  beam  degradation.  Both  approaches  are 
needed,  however  the  turbulence  studies  are  the  more  basic. 

Nearly  all  of  the  propagation  theory  derives  from  the  spatial 
spectrum  of  raicrotemperature ,  as  temperature  turbulence  almost 
entirely  governs  refractive -index  fluctuations  in  the  visible 
to  infrared  wavelengths.  If  the  theoretical  spectrum  were 
incorrect,  the  theories  would  require  modification.  This 
report  will  concern  itself  with  measurements  of  microtemperature 
fluctuations  with  an  emphasis  on  the  low  wavenumber,  input 
range  of  turbulence. 

The  first  indication  of  a  problem  with  the  spatial  spectrum 
was  in  measurements [1]  of  the  phase  difference  power  spectrum 
of  a  spherical  wave  at  10. 6u  wavelength.  The  data,  taken  at 
the  Rome  Air  Development  Center,  demonstrated  more  power  in  the 
low  frequencies  than  was  predicted  by  theory[2].  Even  wnen  the 
theory  is  modified  to  account  for  a  fluctuating  wind  direction, 
the  agreement  is  still  not  good.  (See  figure  1.)  This  led  us 
to  suspect  the  bases  for  the  theory,  since  we  were  confident  in 


Figure  1 •  Experimental  and  theoretical  phase  difference  power 
spectrum  measured  by  Huber  and  TJrt*[l],  Frequency  f  in  this 
and  subsequent  plots  is  in  Hz. 


the  measurement  and  the  processing.  The  bases  are:  the  Rytov 
approximation  in  the  solution  of  the  wave  equation,  the  spatial 
refractive-index  spectrum  and  the  frozen-flow  hypothesis.  The 
Rytov  approximation  should  hold  equally  well  in  the  low  and  high 
frequencies  and  furthermore  is  not  suspect  in  the  calculations 
of  phase  quantities.  Since  the  microtemperature  closely  governs 
the  refractive-index  fluctuations,  we  conducted  a  series  of 
experiments  to  verify  the  von  Karman  spatial  spectrum  and  the 
frozen-flow  hypothesis  for  low  frequencies. 

The  input  range  of  the  spatial  spectrum  is  governed  by  the 
mechanisms  of  turbulence  generation  and  a  beginning  cascade 
process.  Thus  local  terrain  inhomogeneities  are  expected  to 
influence  the  results.  The  aforementioned  experiments  were  run 
under  nearly  ideal  conditions  at  two  sites:  the  Verona.  New 
York  test  site  of  the  Rome  Air  Development  Center  (RADC), 
Griffiss  APB  and  the  laser  test  range  of  Air  Force  Weapons 
Laboratory  (AFWL),  Kirtland  AFB,  New  Mexico.  By  ideal  we  mean 
moderately  high  winds  at  raid-day,  clear  skies  and  high  turbu¬ 
lence  strengths.  Since  the  spatial  structure  function  of 
temperature  is  directly  related,  through  a  Fourier  transform, 
to  the  spatial  spectrum  and  does  not  depend  on  the  frozen-flow 
hypothesis,  we  will  use  the  structure  function  as  a  check  on  the 
validity  of  the  von  Karman  spectrum.  We  will  see  that  the  von 
Karman  spectrum  is  not  the  best  model  for  the  data  and  so  an 


empirical  spectrum  is  developed.  The  model  is  then  properly 
integrated  to  give  the  temporal  spectrum  which  in  turn  is 
compared  with  data.  At  that  point  the  residual  error  can  only 
be  attributed  to  the  frozen-flow  hypothesis.  We  will  see  that 
the  hypothesis  works  remarkably  well  in  the  input  range 
frequencies.  Although  a  specific  spectrum  is  proposed,  the 
mathematics  will  be  kept  sufficiently  general  so  that  it  can 
be  extended  to  analyses  of  non-Kolmogorov  turbulence  or  to 
different  wavenumber  dependencies. 
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II.  THEORY 


Classical  Spectra 

We  begin  by  briefly  reviewing  the  history  of  the  von 
Karman  spectrum  and  the  frozen  flow  hypothesis.  The  most  basic 
form  of  the  temperature  spatial  spectrum  was  developed  in  three 
papers  by  Kolmcgorov[3]  and  via  an  independent  method  by 
Obukhov[i|j  and  applies  only  to  the  inertial  subrange: 

$T  (k)  *  k  “1,/S  A01  «  k  «  A0~\  (1) 

The  scale  lengths  \Q  and  AQ  define  the  limits  of  the  inertial 
subrange  and  are  defined  as  the  internal  and  external  scales 
by  Panchev[5].  Tatarskii[6]  uses  the  symbols  S.Q  and  LQ 
respectively  as  inner  and  outer  scales.  The"  inner  and  internal 
scales  are  identical  but  there  is  a  simple  multiplicative 
relation  between  the  external  and  outer  scales.  The  two 
lengths  will  be  explicitly  defined  after  some  further  develop¬ 
ment. 


Von  Karman  suggested[7»8»9]a  form  for  the  velocity 

-l 

spectrum  for  k  <<  \Q  which  was  later  revised  to  the  temperature 
spectrum[5]  (Panchev’s  sections  i|.5»  6.2  and  6.7): 

ftp  (*0  18  0.0330  Cip2  [x2  +  (1.071/Ao)a  ]  ",lA  (2) 
where  the  numerical  coefficients  are  0.0330  *  5  and 

IBirfi737 
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1.071  *  /  2/ i  \*/2 .  The  introduction  of  the  constant  0.0330 

2  \3rT575T) 

seems  to  be  the  idea  of  Tatarskii[6]  to  reduce  the 
coefficient  to  a  single  strength  of  turbulence  parameter  C^,. 

The  temperature  structure  constant  is  defined  from  the 
inertial  subrange  form  of  the  structure  function: 

Dt  (r)  =  CT2r  2/3  A0«  r  <<A0*  (3) 

The  introduction  of  the  1.071  coupled  with  A0  is  attributable 
to  Reinhardt  and  Collins  [2]  to  make  the  definition  of  AQ  more 
meaningful.  In  fact  the  outer  scale  L0  is  related  to  A0  by 

I>o  *  Ao/I.071  (4) 

in  order  to  make  the  outer  scale  definition  agree  with  the 
generally  accepted  cne  given  by  Strohbehn[10] . 

To  account  for  the  internal  scale  effect  at  the  dissipative 
range,  Golitsin[ll]  (reviewed  by  Panchev[5])  in  I960  suggested 
a  Gaussian  roll-off 

I  1  /  J  2 

*T  (O  ~  K~  exp[-(ic/icM  )  ]  (5) 

where  icH*  5.92/AQ  In  our  current  language  is  the  internal  scale 
wavenumber.  Strohbehn  appears  to  be  the  first  to  make  the 
obvious  connection  between  (2)  and  (5)  to  get 

*_  (k)  *  0.0330  C*  (k*  +  L”2)~xl/‘exp[-(ic/ic  )*]  (6) 

T  TO  M 

Equation  (6)  is  what  is  referred  to  in  the  laser  propagation 

literature  as  the  modified  von  Karman  spectrum.  Since  for  the 

most  part  this  report  is  only  concerned  with  input  and  inertial 
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range  effects,  we  will  concentrate  on  the  valdity  of  the  von 
Karman  spectrun(2) . 


We  are  now  in  a  position  to  define  internal  and  external 
scales  \Q  and  AQ.  The  internal  scale  is  the  separation  r 
at  the  intersection  of  Bq,  (r<<  \Q  )  and  (XQ  «r  <<  AQ). 

This  is  consistent  with  Tatarskii[6]  section  1.13.  At  this 
point  Tatarskii  also  clarifies  the  definition  of  A0  with 
respect  to  the  Kolmogorov  microscales.  The  external  scale  is 
never  clearly  defined  by  Tatarskii,  but  in  essence  Panchev[5] 
uses  the  precise  definition  that  AQ  is  the  separation  r 
at  the  intersection  of  DT  (r<<  AQ)  and  DT  (AQ«  r)* 

Taylor’s  hypothesis [12]  was  the  frozen-flow  assumption 
that  the  turbulence  is  simply  being  transported  by  a  fixed 
point  and  is  not  changing.  It  says  that  the  spatial  and 
temporal  spectra,  after  the  proper  scaling  with  the  transport 
velocity  v ,  are  identical.  Panofsky,  Cramer  and  Rao[13]  state 
the  hypothesis  is  "valid  close  to  the  ground"  for  lag 
distances  up  to  90  m  provided  o^/v  ^  ^  where  a*  is  the 

wind  speed  variance.  This  translates  to  equivalency  of 
spectra  for  temporal  frequencies  f  ^  v/570.  We  do  not  know 
what  is  meant  by  "valid"  or  "close  to  the  ground."  Furthermore 
the  oy/v  criterion  does  not  seem  adequately  restrictive.  We  do 
expect  the  hypothesis  to  hold  in  the  inertial  subrange  for 
fairly  homogeneous  terrain. 
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Theoretical  Relations  Based  on  von  Karman's  Spectrum 

The  first  derivations  are  for  the  statistical  description 
of  the  temperature  fluctuations  and  are  based  on  the  von 
Karman  spectrum(2).  No  internal  scale  effects  are  included. 

The  calculations  and  representative  plots  are  for  the  spatial 
structure  function  and  these  temporal  spectra:  temperature  at 
a  single  probe  and  temperature  difference  at  a  probe  pair.  The 
cross-spectrum  and  coherency  are  then  written  in  terms  of  the 
single-T  and  delta-T  spectra. 

The  temperature  structure  function  is  defined  as 

Dt(p)  «  <(T(*  )-T(?<  ))*>  (7) 

12 

where  the  brackets  Indicate  an  ensemble  average  and  r»|  — 9*  1 1 . 

Homogeneity  and  isotropy  are  assumed  from  the  start.  From 
Tatarskii[6] ,  section  1.13*  the  structure  function  is  written 
in  terms  of  the  three-dimensional  spectrum: 

DT(r)  »  8*/“  [l-sinc(»cr)]  ♦,.,(*}£*£>:.  (8) 

After  integration  of  (2)  in  (8) 

DT(r)«  C*1rl/,[1.0iJ68(r/Lo)‘2/j-0.62029(r/Lo)"l/,K_i/j(r/Lo)] 

(9) 


where  the  numerical  values  are 
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1.0468  =  2A ~  /(3T(5/6) )  and 
0.62029  *  25/J/if  /(3r(5/6)r(l/3) ) . 

The  two  asymptotes  which  are  important  are 

2T<r<<Ao>  -  Cj  r1'*  (10) 

and  Pj,(A  «r)  «  0^  A0’^*.  (11) 

Obviously  the  intersection  of  Di^C*’'<h0)  and  b^.{A0«r)  is  at 

i 

r*A0.  it  was  the  prior  knowledge  of  this  intersection 
that  caused  the  insertion  of  1.071  in  the  spectrum(2).  To 
demonstrate  the  shape  of  the  structure  function,  DT(r)/DT(«) 
is  plotted  versus  r/LQ  in  figure  2.  This  normalized  structure 
function  is  identical  to  1-C(r)  where  C(r )  is  the  autocorrela¬ 
tion  coefficient  of  temperature. 

Since  authors  use  different  definitions  of  spectra — some 
one-sided,  some  bilateral,  some  with  an  extra  2ir  factor — we 
will  start  by  writing  the  spectrum  to  be  used  here  in  terms 
of  the  temporal  covariance  R(t)  in  order  to  define  F(f): 

R(t)  »  /"  F(f )  cos  (2xfx)  df.  (12) 

Tatarskii[6]  in  section  1.20  uses  a  two-sided  spectrum  W(«) 
which  is  related  to  F  by 

F(f )  -  4*  W(2wf).  (13) 

In  that  same  section  Tatarskii  relates  the  temporal  and  spatial 
spectra: 
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*>T(f)  =  8^* v-1  s;.rf/v  *T(<)  KdK. 

Carrying  out  the  straightforward  mathematics  we  find 

PT(f)  =  0.07308  cf  v2/3[f2+f?3~s/e  (15) 

where  f.  *  v _  is  a  break  frequency  and  the  numerical 

2irL 

o 

coefficient  is  0.07308  =  2/C(2ir)a/33r(l/3)].  Figure  3  is  a 
plot  of  a  normalized  FT(f)  versus 

‘fhe  temperature  -  difference  spectrum  F^  is  now  calculated 
in  a  manner  somewhat  paralleling  the  phase-difference 
calculations  of  Reinhardt  and  Collins[2]  and  of  Clifford [l1! ] . 
Inverting  equation(12) : 

FAt  (f;ft)«2/^  dt  cos(2uTt )  R^T(r;p)  (16) 

where  p  Is  the  vector  separation  of  the  two  probes.  After 
using  some  geometry  and  the  frozen- flow  hypothesis,  we  obtain 

RA?(T*^)  *  -DT(v*>  +  7  «¥*•*>  *■  i  Dt(p-vt).  (17) 

If  we  assume  the  temperature  field  Is  homogeneous  and 
isotropic  and  that  the  angle  between  p  and  v  is  0,  then 

RAT(T;r)““DT  r  dt  (/(r+vTcose)*  +(vtsin0)i)  + 

-*•  DT(/(r“VTCose)1+(VTsine)1) ,  (18) 

where  r*l|p||«  Inserting  equation  (8)  into  (18)  into  (16) 
and  interchanging  the  order  of  integration: 
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FAT( f ;r)  =  8ir  d<  4>T(tc)<2  £  drcos  2irfT 
[2sinc(fcvx)-sinc(ic»^t’+VTcos0)2+{vTsin0)2 ) 


Having  done  the  T-integration,  we  obtain 

F.m(f  ;r)»l6ir2v  1  /*  d»c4>rr(»c)ic 
Ai  airf/v  1 

[1-cos  ( 2irfrv~ 1  cos0)  JQ(  rsin0/ic2-(2irf/ 


(20) 


After  substitution  of  the  von  Karman  spectrum(2)  and  after 
replacing  u2+l«  (§^f : 

FAT^f»r^*  9ffl737  ctv  X^2irf/V^ 

r0  du  uCu*+l+(f  /f)*r‘*/*[l-cos(fS2tA.)  jo(piS_i 


where 

f  «  _5L_ 

i  2*L0 

and  f  *  . 

*  2irr 


(21) 


(22) 


(23) 


are  the  characteristic  frequencies.  After  final  integration 


of  (21): 


FaT( *>3.1272  CTr5Av"1Sle-;/2  [i'cos  (1'6(|J)] 

where  i.g(B)  -t(f|Arai7^  Kf(6)-  (25) 

82  .  [l.(£)J](pi"<>)  !  .  '<»> 


(26) 


G(0)»0 


(27) 


and  the  numerical  coefficient  is  3«1272»2^ir/(9r(l/3)). 


^  S’} 


The  function  K  (0)  is  the  modified  Bossel  function  of  the 
v 

second  kind,  vth  order[15]*  Figure  ll  Is  a  plot  of  the 
normalised  F^T(f)  versus  f/f^  for  f2/x‘i=Lo/r=10 .  The 
resonances  in  the  high  frequencies  f>f2  occur  for  0  near  0 
or  ir  since  similar  turbulence  is  seen  at  both  sensors  when  the 
wind  is  blowing  along  the  sensors.  Since  the  time  lag  of 
such  a  flow  is  rv”1cos0,  the  minimums  in  the  spectrum  occur 
at  fjj»Nvr-1  sec-0,  where  N  is  an  integer.  In  the  real 
atmosphere,  however,  the  wind  speed  and  direction  are  not 
the  fixed  values  required  by  fN.  The  small  variations  that 
normally  occur  in  v  and  0  are  sufficient  to  wash  out  the 

theoretical  resonances.  For  example  when  0*0  and  f>>f2,  the 
1-cos (  )£l.  It  is  interesting  to  consider  the  two  regions 

which  have  simple  power-law  dependencies  for  0*ir/2: 
FAT(f«fJf2)*5.829  C*  r5/sv-»  (28) 

and 

FAT(fif2«f)  *  0.14616  v2/3C2f“5/3.  (29) 

Hence  the  difference  spectrum  approaches  the  constant  given 
by  (28)  as  f-*o,  and  has  the  high  frequency  dependence  (29) 
which  *s  exactly  twice  the  inertial  range  dependence  of  FT 
given  by  (15)  for  f»f^: 

PAT(fi»f2<<f)  *  2FT(f1«f).  (30) 
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In  a  similar  fashion  there  are  three  regions  when  0=0  or  ir: 

* 

FAT(f<<f\’  V'3*-1272  CT  Lo5/3v-‘  T  (3D 

F.m(f  «f«f  )=3.1272  4  r^v"1  ~  (f/f  )l/3,  (32) 

and  FAT(fi«f)»3.1272  C|  rs/sv"1  2(f/f 2)“s/3sin2 (33) 

Hence  the  power  law  dependencies  are  f 2 ,  f1^3  and  f~5/^3  in  the 

three  ranges.  If  the  wind  velocity  is  sufficiently  fluctuating 

i 

that  sin2(  )  be  replaced  by  its  average J2  ,  then 


F  T(f^,  f2<<f)^°.  1^6.16  v2/jC2  f"s/3. 


(3*0 


which  is  the  same  as  equation  (29). 

The  remaining  spectra — the  cross-spectrun  and  the 

coherency — are  easily  calculated  from  the  single-T  and  the 

delta-T  spectra.  The  cross-spectrum  F  is  complex  in 

12 

general  and  has  as  Its  real  part  the  cospectrum  CT  T  and  as 


its  imaginary  part  the  quadrature  Qrj.  >p  : 


1  2 


1  2 


FT  t  (?)  *  CTjT2^  +  i^TiT2^ 


(35) 


1  2 


The  coefficient  of  coherency  is  a  sort  of  ’’spectral  correlation 
function”  and  is  in  general  complex: 

PT  t  (f) 

v  (f)  «  - — - -  •  (36) 

V.  /?T|(»F~TfT 


The  phase  of  the  cross-spectrum  and  the  phase  of  the  coefficient 
of  coherency  are  identical: 


9TiT2(f)sStan~1(QTiT2(f)/CTiT2(f  ^ 


(37) 


The  phase  relates  any  temporal  shift  of  one  signal  T  with 

respect  to  a  second  T  .  The  modulus  of  the  coefficient  of 

2 

coherency,  termed  here  simply  as  the  coherency,  is  a  measure  of 
how  one  signal  relates  to  another  as  a  function  of  frequency.  A 
perfect  relation  is  unity  and  no  similar  trend  at  that  frequency 


is  zero;  thus 


0£  I  kTiT2  I— 


(38) 


If  the  reader  desires  a  further  reference  to  these  quantities, 
he  might  see  Robinson[l6] . 


The  previous  paragraph  was  a  general  treatment  of  the  cross¬ 
spectrum  and  coefficient  of  coherency  for  two  signals  T^  and  T^ 
which  in  the  present  context  are  temperatures  measured  at  two 
probes.  We  will  return  to  these  quantities  in  the  data  sections. 
However,  the  theoretical  individual  spectra  are  assumed  equal. 


so  the  cospectrum  (plotted  in  figure  5)  is 


CTiT2<f>  *  FT(f)-f  FAT(f) 


(39) 


The  phase  is  formed  easily  from  geometrical  considerations  as 
f  4.  cose  rtraosce 

9(f)*  4  U  .  C'Oa) 

I  0  f>f3csc0. 

The  theoretical  coherency  (figure  6)  is  then 


lT,T* 


(f)|«  (l  - 


2FT(f) 


sec@. 


(40b) 


% 


III.  DATA  ACQUISITION  SYSTEM 


The  Recording  Scheme 

The  purpose  of  the  data  acquisition  and  processing  system 
is  to  provide  temporal  and  spatial  statistical  information  on 
the  turbulent  temperature  of  the  atmosphere.  The  overall  system 
is  diagrammed  in  figure  7-  Common  meteorological  observables 
which  are  recorded  include  wind  speed,  wind  direction,  baro¬ 
metric  pressure  and  the  gross  temperature  from  a  dew-point 
system.  A  description  of  all  the  sensors  will  follow.  Nine 
microthermal  probes  are  oriented  in  a  geometrical  progression 
on  a  horizontal  boom  posltionable  two  to  eight  meters  above 
the  ground.  This  makes  a  total  of  thirteen  signals  which  were 
recorded  digitally  on  a  narrow  bandwidth  system  and  on  analog 
tape  for  subsequent  wide  bandwidth  processing.  The  digital 
unit  is  the  Analog-Digital  Data  System  model  100m  Data  Logger, 
with  a  12  bit  analog-to-dlgital  converter.  The  sample  rate 
is  **00  Hz  divided  by  the  number  of  channels  currently  in  use. 

For  the  13  channels,  the  rate  is  about  31  Hz.  This  narrow 
bandwidth  system  is  used  for  obtaining  those  statistical 
quantities  which  are  not  sensitive  to  aliasing:  especially 
probability  densities,  variances,  spatial  cross-covariances 
and  higher-moments.  For  spectral  information  on  the  turbulence, 
all  thirteen  signals  are  recorded  in  parallel  on  analog  tape. 

The  Ampex  FR1200  tape  recorder  has  at  3  3/1*  ips  a  bandwidth  of 
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DC  to  1250  Hz  (-3  db)  and  the  tapes  used  have  bandwidths 
superior  to  that.  Later  two  channels  at  a  time  are  converted 
from  analog-to-digital  on  a  12-bit  PDP-8  computer.  Because 
the  possibility  of  aliasing  will  be  significant  to  the  subse¬ 
quent  Fourier  analysis,  a  low-pass  filter  is  inserted  in  each 
channel  prior  to  conversion.  The  -3db  point  of  the  Krohn- 
Hite  model  3322  filters  is  set  to  the  Nyquist  frequency,  one- 
half  the  sample  rate,  and  the  roll-off  is  that  of  a  four-pole 
Butterworth  filter,  -80db/decade.  Sample  rates  from  400  to 
700  Hz  are  chosen  depending  on  the  expected  noise  level.  Itoe 
measured  Krohn-Hite  power  response  is  plotted  in  figure  8 
for  a  cut-off  frequency  of  200  Hz.  To  estimate  the  response 
at  higher  cut-off  frequencies,  you  merely  slide  the  curve  to  the 
right  until  the  -3db  point  is  at  the  desired  frequency. 

Digital  Processing 

The  digital  tapes  from  the  ADDS  unit  and  from  the  PDP-8 
are  processed  in  two  separate  software  packages  on  a  Honeywell 
635/645  computer.  One  program  computes  the  cross-covariance 
matrix  of  the  thirteen  signals  by  operating  on  the  tapes  from 
the  ADDS  unit.  Matrices  are  generated  for  averaging  times 
typically  ranging  from  10  to  30  minutes.  The  microtemperature 
elements  of  the  matrix  are  then  converted  into  a  plot  of 
l-C(r)  where  C(r)  is  the  spatial  correlation.  Such  a  plot  is 
identical  to  DT(r)/Drr(«) ,  found  theoretically  for  the  von 
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Pigure  8.  Power  responses  of  aliasing  filters  (2  channels) 
and  the  expeoted  response  of  2.3m  sensors. 


Harman  spectrum  in  the  previous  section.  In  the  end  an  average 
l-C(r)  plot  is  generated  for  the  entire  mission.  Overall 
averaging  times  vary  and  will  be  Indicated  in  the  mission 
summaries  yet  to  come.  This  program  also  calculates  Ci| 
according  to  its  definition  (7)  and  calculates  AQ  as  the 
separation  r  where  the  extrapolated  inertial  subrange  of  l-C(r) 
intersects  the  unity  line. 

The  second  program  operates  on  the  tapes  from  the  PDP-8 
computer  and  provides  wide  bandwidth  power  spectral  densities. 

The  six  spectra  were  described  in  the  previous  section:  the 
individual  spectra  of  the  two  temperatures,  the  difference 
spectrum,  the  modulus  and  phase  of  the  cross-spectrum  and  the 
coherency,  The  signals  are  first  transformed  via  a  fast  discrete 
Fourier  algorithm  for  frequencies  from  fny/20ll8  to  fny,  where 
f ny  is  the  Nyquist  frequency.  Namely  *4096  data  points  per 
channel  are  transformed  at  one  time.  The  data  are  then 
digitally  low-pass  filtered  by  a  four-pole  Butterworth  with 
a  cut-off  of  fny/6^.  The  low-frequency  part  of  the  spectrum 
is  then  obtained  by  transforming  the  remaining  signal  which 
has  been  decimated,  keeping  only  one  point  in  sixty-four. 

The  low  and  high  frequency  spectra  undergo  constant  percentage 
bandwidth  averaging  so  that  on  a  logarithmic  frequency  scale, 
the  spectral  estimates  appear  to  have  a  nearly  constant 
separation.  Sixty-four  high  frequency  spectra  are  averaged 


to  correspond,  in  averaging  time,  to  a  single  low-frequency 
spectrum.  Finally  several  such  spectra  are  averaged  to 
represent  an  entire  mission.  The  final  averaging  time  can  be 
computed  as  W2048«6Vfny ,  where  N  is  the  number  of  spectra 
averaged  together.  Since  varying  numbers  of  spectral  estimates 
are  averaged  to  give  a  final  value,  the  confidence  limits  are 
a  function  of  frequency  as  well  as  N.  Following  Blackman  and 
Tukey[17],  we  assume  the  spectral  estimates  follow  a  chi- 
squared  distribution.  Tables  in  Blackman  and  Tukey  or  in 
Abramowitz  and  Stegun[15]  can  be  used  to  generate  upper  and 
lower  confidence  limits.  The  confidence  limits  for  our 
spectra  never  exceeded  twice  the  value  at  the  high  side  or  1/2 
the  value  at  the  low  side  and  the  error  bar  size  decreased 
rapidly  with  increasing  frequency. 

Commercially  Available  Sensors 

The  standard  meteorological  sensors  used  are:  (a)  Packard- 
Bell  wind  speed  indicator  in  a  "staggered  six"  configuration, 
model  52.3  with  a  model  50.1  transmitter;  (b)  Packard-Bell 
wind  direction  indicator  "quick-one"  model  53-1  with  a  model 
50.2  transmitter;  (c)  Sostman  pressure  transducer  model 
20114-28/32HA-l  with  a  model  1230  signal  converter  and  (d) 
Cambridge  Systems  Inc.  dewpoint  system  model  110-S-M.  These 
systems  generate  low-frequency  signals  and  are  not  expected  to 
correlate  with  the  fast-response  temperature  measurements. 
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Mlcrothermal  Sensors 


While  there  are  commercial  systems  available, the  fast- 
response  microthermal  sensors  and  amplifier  circuits  are 
constructed  in  the  laboratory  because  of  significantly  lower 
costs  and  a  better  understanding  of  what  the  instrument  consists 
of.  The  sensors  and  amplifiers  are  based  on  a  design  of  Ochs 
at  N0AAC18].  The  sensors  themselves  are  simply  a  short  piece 
of  Wollaston  process  silver  on  platinum  wire  which  has  been 
etched  with  HNO^  to  give  an  active  element  2mm  long  and  about 
2.5v  in  diameter.  Prior  to  etching,  the  wires  are  connected 
with  cold  solder  to  the  supports.  While  our  supports  are 
specially  constructed,  an  old  light  bulb  base  works  well  and 
gives  a  convenient  way  to  swap  sensors.  The  small  platinum 
wires  sense  temperature  changes  via  resistive  changes.  The 
wire  is  simply  placed  as  one  leg  of  a  Wheatstone  bridge  and  the 
imbalance  voltage  is  a  function  of  resistance.  Since  only 
about  450pa  current  is  drawn  through  the  wires,  they  are 
effectively  cold-wire  probes  as  long  as  the  wind  velocity 
exceeds  about  0.5  m/sec. 

0chs[l8]  uses  sensors  which  arc  0.6p  in  diameter,  and  they 
have  a  rvspor.se  time  of  about  200psec  in  still  air  or  about 
70-80psec  in  winds  of  2  to  5  m/sec.  Measurements  of  the 
response  time  of  a  2.5u  wire  were  made  by  Ochs [19 ]  as  900psec 
in  still  air  and  300psec  in  winds  of  about  5  m/sec. 
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Wyngaard[20]  corroborates  tnese  values  with  measurements  of 
the  response  time  of  a  sensor  used  in  the  RADC  experiments, 
finding  about  800ysec  in  still  air.  The  values  are  for 
"dirty"  sensors,  ones  which  have  been  used  for  a  day  or  so. 
"Clean"  sensors,  freshly  made,  may  have  response  times  at 
least  one-half  the  aforementioned  values.  By  response  time 
we  mean  the  time  it  takes  the  output  voltage  of  the  bridge 
to  go  to  63. 2%  of  its  final  value  after  the  sensor  sees  a 
step  change  in  temperature-.  Since  the  response  is  exponential > 
the  sensor  acts  as  its  own  low  pass  filter  with  a  power  response 
of 

[1+  (2irfT)2]-1  (41) 

where  T  is  the  response  time.  This  is  plotted  in  figure  8 
for  T=300  and  900psec.  The  length  of  the  sensor  determines 
the  minimum  spatial  scale  that  the  probe  can  respond  to.  As 
our  probes  are  2mm  in  length  we  can  expect  accurate  measurements 
of  statial  structure  functions  for  scales  down  to  about  6mm  and 
spectral  measurements  to  about  f»v/(2ir0.006).  Unfortunately 
precise  observations  of  the  internal  scale  will  not  be  possible 
since  we  expect  XQ^5mm. 

To  summarize  the  effects  of  various  low-pass  filters  at 
the  high  frequency  end  of  the  spectra  we  have:  (a)  the  cutoff 
of  the  aliasing  filter  which  is  so  sharp  as  to  not  have  an 
effect  below  0.7fny;  (b)  the  300ysec  response  time  of  the 
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sensor  which  gives  a  -3db  point  at  530Hz  and  has  essentially 
no  effect  below  about  170Hz;  (c)  the  physical  cut-off  of  the 
turbulence  in  the  dissipative  range  which  is  approximately 
Gaussian  with  a  e-1  point  of  f=v/(2irA0)  and  (d)  the  length  of 
the  sensor  which  has  an  unknown  effect  in  averaging  scale 
sizes  on  the  order  of  2mm  and  smaller  but  probably  has  a  cut¬ 
off  frequency  around  f=v/(27i0.002) .  If  all  this  does  not 
sufficiently  muddle  the  high  end  of  the  spectrum  for  the 
reader,  he  may  also  correctly  suspect  that  at  this  point  the 
spectrum  is  about  to  go  Into  the  noise.  At  the  low  end  of  the 
spectrum  there  is  a  single  high-pass  filter  to  remove  the 
effects  of  diurnal  temperature  changes  which  are  not  low 
frequency  turbulence  effects.  The  lOOsec  time  constant  of 
that  simple  RC  filter  gives  a  -3db  point  at  1.6xlO_JHz. 

The  bridge/amplifier  schematic  diagram  Is  given  as  figure 
9  and  the  component  list  Is  table  1.  After  the  sensors  have 
been  constructed,  their  resistance  is  measured  and  scaled  to 
0°C.  The  Calender-Van  Dusen  scaling  equation  as  given  by 
Rosemount  Bulletin  134l[21]  applies  to  pure  annealed,  strain- 
free  platinum  wire  in  the  temperature  range  -183£T<630°C: 

gg&P-  .  1  +  o[T+ {(l-T/100) (T/100)  +  6(1-T/100) (T/100) s] 

KsVO  ) 

(42) 


TABLE  1:  COMPONENT  LIST  FOR  BRIDGE/AMP LIFIEP  CIRCUIT 


CODE 

QUANTITY 

VALUE 

A1 

1 

- 

A2 

1 

- 

Cl 

1 

5CGpfd 

C2 

2 

200yfd 

C3 

1 

390yyfd 

C4 

1 

220yyfd 

C5 

1 

• 

lQyyfd 

Ra 

C 

30.  IK 

Rb 

X 

2.7^ 

R 

l 

3 

1  Meg 

R 

2 

3 

1  K 

R 

3 

1 

10K 

RP 

1 

200fi 

2 

- 

- 

2 

- 

rs(t) 

1  spool 

- 

- 

1 

- 

1 

ITEM 

AD208H  Op.  Amp. 

AD504LH  7235H  Op.  Amp. 

15  WVDC  Mallory  MTV  500  DW15  Capacitor 

25  WVDC  Sprague  TE-1213,  Capacitor 
30D207G025DH2 

±5%  Elmenco  Mica  7FA391G033  Capacitor 

± 5 %  Elmenco  Mica  7FC221G03  Capacitor 

±5?  Elmenco  Mica  6CD100K03  Capacitor 

±1?  TI  Resistor  MC65C 

±1S  TI  Resistor  MCfi5C 

±1?  TI  Resistor  MC65C 

±1%  TI  Resistor  MC65C 

±1%  TI  Resistor  MC65C 

Amphenol  Potentiometer  99^SL201 

Amphenol  Connector  7  Pin  #197 

IC  Socket  Cinch-Jones  8-ICS 

Wire  0.1  mil  Diam  Platinum  Silver 
Coated,  57^7«/Ft,  20°C,  #408, 
Sigmund  Cohn  Corp,  Mt.  Vernon,  NY 

Power  Supply  Zeltex  ZM1550 

Subminiature  Switch,  Switchcraft 
JMT-121 
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where  <?  =  0.00393 
6  =  1.49 

B  *  0.110  for  T<0  and  8*0.00  for  T>0. 

The  sensor  Indicated  by  R-  is  installed  as  one  leg  of  the 
bridge  and  the  potentiometer  plp  is  adjusted  such  that  the 
output  voltage  is  nearly  zero.  The  gain  equation  of  the  entire 
bridge/  amplifier  system  is  calculated  for  better  than  155  error 
as 


Eout<T> 


T-T-, 


E  G  R_ (0) 

a  (1+6/100) 

^T~ 


(43) 


where  *  temperature  (°C)  when  bridge  is  nulled, 

E  *  bridge  bias  voltage  (-15vt) 

Rb  «  bias  resistance  (2.74K) 

Ra  *  resistances  in  two  legs  of  bridge (30. IK) 
and  G  *  voltage  gain  of  amplifier  (-l^xlO'’). 

Substituting  in  these  values: 


WT>  V°> 

M 

T-T^  54.1 


(44) 


where  Rg(0)  is  in  ohms  and  is  ln  Vt/°C.  Hence  a 

sensor  resistance  of  54.10  will  yield  a  lvt/°C  gain  of  the 
system.  To  convert  the  actual  output  voltage  to  temperature 
changes,  multiply  Eout  by  54.1/RS(0). 


t 
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Amplifier  Drift  Test 


Since  we  will  be  presenting  spectral  densities  into  the 

low  frequencies,  below  1Hz,  we  are  concerned  with  electronic 

component  drift  in  the  bridge/amplifier  circuit.  To 

empirically  check  the  noise  level  of  the  system,  we  placed 

a  micro thermal  probe  in  a  bath  of  2  qt.  oil.  Except  for  the 

oil  bath,  the  system  was  run  normally.  The  oil  temperature 

had  been  allowed  to  stabilize  to  the  ground  temperature  before 

the  mission.  A  second  sensor  was  placed  in  air  approximately 

two  meters  above  the  sensor  in  oil  in  order  to  provide  the 

turbulence  statistics .  The  spectral  data  obtained  from  these 

two  sensors  on  28  Jun  73  are  presented  in  figure  10.  Total 

mission  time  was  20  min.  Even  for  this  low  turbulence  day 

when  C*  *  0.004  °C2m“2/s,  the  noise  level  is  at  least  20db 
T 

below  the  signal  spectrum,  except  at  the  highest  frequencies. 
The  low  frequency  power  is  attributed  to  amplifier  drift  and 
oil  temperature  variations.  The  series  of  spikes  around  2  to 
5  Hz  came  from  the  recorder  heads.  This  problem  has  since 
been  corrected.  The  remainder  of  the  power  is  the  composite 
noise  of  the  bridge, amplifier,  coaxial  cable,  tape  recorder 
amplifiers  and  analog  tape.  The  conclusion  is  that  drift  and 
other  r.oise  sources  do  not  significantly  influence  the  data. 


IV.  DATA  COMPARED  WITH  VON  KARMAN'S  MODEL 


We  now  investigate  the  data  of  a  typical  data  collecting 
mission  at  Grif'fiss  APB  to  show  the  regions  of  agreement  and 
disagreement  with  the  theoretical  curves  of  section  II  which 
were  based  on  the  von  Karman  spectrum.  The  Verona  test  site 
near  Griffiss  is  where  the  predominance  of  data  is  taken. 

Verona  is  in  the  center  of  the  Mohawk  Valley,  essentially  flat 
terrain.  The  local  terr'ain  is  best  shown  in  a  site  sketch, 
figure  11a.  Prevailing  winds  are  typically  around  285°  or 
about  135°-  Mission  210873(21  Aug  73)  was  run  at  mid-day 
during  the  most  unstable  conditions  of  turbulence  with  a 
clear  sky  and  moderate  winds.  The  conditions  of  the  mission 
are  best  displayed  in  the  mission  summary,  figure  lib.  This 
will  be  the  standard  format  for  other  missions  to  be  presented. 
The  four  plots  are  mean  wind  speed  v  and  its  standard 
deviation  ov,  mean  wind  direction  3  and  its  standard  deviation 
o q  with  respect  to  the  line  of  sensors,  the  external  scale 
A0,  and  the  temperature  structure  constant  squared  C^>,  all 
versus  time  during  the  mission.  The  standard  deviations 
are  shown  as  X's  on  the  same  plots  as  the  mean  values. 

The  parameters  needed  to  develop  the  theoretical  curves 
are  also  presented  in  figure  lib.  These  theoretical  curves 
and  the  spectral  and  spatial  data  are  presented  In  figures  12a 
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Figure  11a.  Site  sketch,  Verona  test  site  of  Griff iss  AFB  KY. 
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through  12f.  In  figure  12a  the  averaged  normalized  structure 
function  is  plotted  versus  separation  r  for  0.04<r<6m.  In  the 
large  separations,  when  r^0.6m,  the  data  noticably  deviate 
from  the  theoretical  curve.  This  means  there  is  a  better 
correlation  between  the  temperatures  on  probes  in  the  vicinity 
of  the  external  scale  than  is  predicted.  When  translated  into 
optical  propagation  terms  this  means  a  wider  coherence  function 
than  previously  predicted.  The  data  further  indicate  that  a 
modified  curve  is  needed  which  would  depart  from  the  inertial 
subrange  at  a  smaller  separation  and  have  a  more  gradual  knee 
in  approaching  the  unity  asymptote.  In  the  least  the  data 
demonstrate  the  turbulence  was  Kolmogorov,  with  an  r2^3 
dependence  in  the  inertial  range. 

The  spectra  in  figures  12b  through  12f  are  for  two  sensors 
selected  such  that  their  separation,  r  =  0.289m,  is  the 
inertial  range.  Figure  12b  is  a  plot  of  the  power  spectra  of 
the  individual  temperature  signals.  Apparently  the  two  sensors 
statistically  saw  essentially  the  same  fluctuations.  The 
theoretical  curve  there  shows  a  discrepancy  in  the  low 
frequencies  for  f<?f^.  This  could  have  been  predicted  by  figure 
12a  since  large  separations  correspond  to  low  frequencies. 
Agreement  in  the  inertial  subrange  is  excellent  and  adequately 
shows  the  Kolmogorov  power  law,  f"5/3.  In  the  frequencies 
10£f<100Hz  there  is  slightly  more  power  than  an  inertial 
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Figure  12b*  Normalized  temperature  spectral  data  and  theoretical 
curve  based  on  von  Karmen's  spectrum*  Individual  spectra  were 
normalized  prior  to  averaging  for  this  final  curve* 


Figure  12c.  Formalized  temperature  difference  spectral  data  and 
theoretical  curve  based  on  von  Kansan*  a  spectrum. 


subrange  allows  '’or.  This  can  only  be  in  the  turbulence 
itself  and  is  presently  inexplicitable  except  to  the  extent 
that  this  region  of  the  spectrum  is  influenced  by  the  internal 
scale . 

The  temperature  difference  spectrum  for  mission  210873  is 
in  figure  12c  and  there  the  excess  low  frequency  power  is  no 
longer  attributable  to  an  error  In  the  von  Karman  snatial 
spectrum.  In  fact  a  difference  measure  for  ©£tt/2  should  be  a 
fairly  smooth  flat  curve  below  f=f2.  The  sensor  pair  is  not 
sensitive  to  changes  in  the  distribution  of  turbulent  eddy 
sizes  for  scales  significantly  longer  than  the  separation. 
Rather,  the  excess  power  seen  in  the  data  is  attributable  to  an 
inappropriateness  of  the  frozen  flow  hypothesis.  The  value  of 
av/v  for  this  mission  is  0.27  ,  so  it  is  not  surprising  that 
the  hypothesis  is  not  working  well  for  scale  sizes  longer  than 
18m  in  the  difference  spectrum. 

Figure  12d-e  is  the  cross-spectrum,  obtained  directly  from 
the  raw  data  and  not  from  the  difference  spectrum.  For  those 
frequencies  below  f2,  where  the  difference  spectrum  is 
significantly  below  the  individual  spectra,  the  cross-spectrum 
approximately  equals  the  individual  spectra.  For  the 
frequencies  above  f2  the  turbulent  eddies  appear  independent 
at  the  two  sensors  and  so  the  cross-spectrum  goes  into  the 
the  noise  quickly.  Theoretically  the  Dhase  is  given  by  (JJOa). 
The  spread  in  the  data  above  f%3Hz  is  due  to  a  lack  of 
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averaging  a  sufficient  amount  of  data. 

Finally,  figure  12f  is  the  coherency.  The  coherency  and  its 
associated  phase  angle  in  the  previous  graph  show  how  much 
aliKe  the  two  signals  are  and  are  excellent  tools  in  inuring 
the  processing  is  done  correctly  in  the  high  frequencies.  For 
scale  sizes  longer  than  the  separation  the  coherency  is  near 
unity,  but  for  scale  sizes  increasingly  smaller  than  the 
separation,  the  coherency  decreases  until  it  reaches  the  noise. 
It  is  interesting  that  the  noise  above  10Hz  on  the  two  channels 
has  a  relatively  high  coherency.  Thus  the  noise  is  coming 
from  the  same  source,  possibly  the  tape  recorder  heads  or  the 
analog  tape. 

In  summary  there  is  obviously  some  inaccuracy  in  the 
von  Karman  spatial  spectrum  for  wave  numbers  in  the  input 
range.  The  errors  were  translated  into  discrepancies  between 
theory  and  data  in  the  structure  function  for  r^A0and  in  the 
spectra  for  f<f^.  The  following  section  will  develop  an 
empirical  model  for  the  spatial  spectrum  which  better  matches 
the  data.  Since  one  mission  is  not  sufficient  data  to 
establish  such  a  model,  we  will  then  show  how  the  new  theory 
compares  with  the  data  of  three  other  missions  at  two 
radically  different  sites  and  at  three  altitudes. 


V.  EMPIRICAL  MODELING 

At  present  we  do  not  know  how  much  of  the  excess  low 
frequency  power  is  attributable  to  the  use  of  von  Harman' s 
spectrum  and  how  much  to  the  invalidity  of  Taylor's  hypothesis. 
Since  the  structure  function  is  purely  a  spatial  measurement 
not  depending  on  frozen  flow  and  since  it  is  related  directly 
to  the  spatial  spectrum  through  (8),  it  will  be  the  basis  for 
the  development  of  a  new  model  for  the  spatial  spectrum.  The 
model  will  be  then  treated  as  was  the  von  Karman  form  in  that 
the  temporal  spectra  will  be  derived  and  compared  with  the 
data. 

General  Form  of  the  Model 

We  assume  the  form  of  the  spatial  spectrum  is 

*T(0  *  c[(<L0)q  +  (<L0)S]"P,  (45) 

where  once  again  th*»  internal  scale  effects  are  ignored.  The 
coefficient  c  will  soon  related  to  C^i  and  to  LQ.  This  eouation 
reduces  to  the  von  Karman  form  when  q*2,  s=0 ,  p=ll/6  and 
c=0.0330  Lq1/,j.  The  condition  for  finite  power  is  that 

Dt(®),  which  is  twice  the  temperature  variance,  be  finite: 

Dfp(®'  *  8ir/~  ♦T(ic)K2dic«».  (l}6) 

This  is  satisfied  if  sp<3<qp,  which  in  turn  implies  the 
power  law  dependence  -sp  in  the  low  wavenumbers  is  always 
greater  than  the  high  frequency  dependence  -qp.  This  insures 
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a  knee  in  the  curve  in  the  vicinity  of  k'v-Lq”1.  After 
substituting  (45)  into  (46)  we  get 


Dip  (® )  =  8-itc  1 


fci2) 


We  will  not  attempt  to  adjust  c  so  that  D^t®)  is  not  a 
function  of  s,p  and  q.  Rather  c  is  set  bv  the  inertial 
subrange  form  of  D^,(r).  Since  Dip(r<<LQ)  does  not  depend  on  the 
low  wavenumber  shape  of  and  since  the  total  power  is 
finite,  we  can  calculate  DT  (r<<Lc)  by  temporarily  assuming  a 
form  much  simpler  than  (45): 


4>t(ic)  *  c[kL0+1]-c*P  •  (48) 

After  substituting  (48)  into  (8)  we  obtain  a  confluent 
hypergeometric  function  which  is  then  expanded  about  r^O. 
Keeping  only  the  lowest  order  term: 

DT(r<<LQ)=c  rP9"JL0”pci8ir2 . 

cscy?-pq.U  sln  [ir  ( pq-2)/2 1  (49) 

f(pq-l) 

Even  though  this  is  a  direct  result  of  using  (48),  this 
equation  applies  for  all  pq>3*  When  pa  is  sin  integer, 
L'Hospital's  rule  must  be  used  to  evaluate  (49).  In  keeping 
with  the  accepted  difinition  of  in  (10),  we  let  CiJ2  be  the 
constant  of  proportionality  in  the  inertial  subrange: 


DT<r«L0)«C*2rp(3-J. 


(50) 
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(51) 


A  comparison  of  ( ^49 )  with  (50)  reveals 
C*2 

c—  L?q  r(pa-l)sinO(2-pq)  3csc[ir(pq-2)/2]  . 

8u2  0 

Another  important  quantity  which  is  derived  from  the 
general  spectrum  is  the  external  scale  length  AQ.  In  the  case 
of  a  modified  von  Karman  spectrum  Ao=1.071Lo.  Once  again,  the 
external  scale  is  the  value  of  r  where  D<j>(r<<L0)/D.j>(a’)=l.  Vie 
evaluate  Ac  from  equations  (47)  and  (49): 


iA0-.pq-3  r(^!r(^)rtP^> 

K> 


(q-s)irr(p) 

sin[t(2-pa)]csc[if(pq-2)/2].  (52) 

As  for  the  general  temporal  spectrum,  we  can  derive  the 
low  and  high  frequency  asymptotes.  By  substituting  (45)  into 
(.14)  we  find  for  op >2  ,  0<sp<2  and  f<<f^  : 


vL2  r(£^£}iiS2z2\ 

Lo  F<ti(  f )  •*■  1  l  q-s  I  \q-s  / 

T(p) 


,  ^  ,2-sp 

SrT ) 


bTT2  C 


d-S 


2 -so 


(53) 


and  for  2<qp 


and  f>>f- 


vLo  Fm(f)  _ 

8itJc  ‘  P9-2  [rl' 


IU 


2-pq 


(54) 


The  requirement  that  sp<2  only  applies  to  equation  (53)- 
Certainly  the  power  is  still  finite  for  2<sp<3,  although  now 


iim  Fm(f A  general  form  for  the  difference  spectrum 
f-*0 


F.m(f)  is  quite  difficult  to  obtain  analytically,  however  we 


know  that 


FiT(fI’  f2<<f)  =  2FT(fi<<f)‘  (55) 

Introducti on  of  a  Kolmogorov  Iner t i a I  Range 

For  the  remainder  of  this  paper  the  power  law  dependence 
in  the  inertial  subrange  will  be  Kolmogorov.  That  is, 
qp=ll/3-  Hence 

0=0.0330  C£  L0,l/2.  (56) 

We  are  more  concerned  with  the  power  spectra  at  lev;  frequencies 
and  the  structure  function  at  large  separations  than  with 
small  differences  in  the  inertial  range.  The  normalised 
spatial  spectra-. -^(ic)  versus  is  plotted  in  figures  13a, 

13b  and  13c  for  q=l,  q=2  and  a=4  respectively.  The  curves 
demonstrate  that  a  adjusts  the  rate  of  transition  from  the 
input  range  to  the  inert ial  range.  The  higher  q's  give  a  more 
abrupt  transition.  The  value  of  s  used  changes  the  power  law, 

-sp,  in  the  low  wavenumbers. 

Numerical  integration  is  required  to  generate  DT(r)  for 

any  q  an  1  sp,  the  two  governing  parameters.  Since  Dr^®)  is  U 

\  7 

given  by  equation  (47),  we  need  only  to  plot  DT(r)/DT(®)  • 

versus  r/L0,  and  this  is  done  in  figures  l4a-c  with  the  same 
parameters  as  in  figures  13a-c.  The  best  fit  of  these  curves 
to  the  data  of  mission  210873  is  found  by  overlaying  that  data 
and  by  sliding  the  experimental  curve  along  the  abscissa.  ’ 
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sp=1/2 


It  is  rather  apparent  that  the  q=l  curves  have  too  gradual 
a  transition  between  the  inertial  and  the  input  ranges  and  that 
the  q=1l  curves  are  too  abrupt.  The  parameters  q=2  and  sp=ll/6 
seem  to  give  the  optimum  match.  In  figure  15  the  structure 
function  data  of  mission  210873  is  overlayed  on  that  optimum 
curve.  Least  squares  fitting  to  obtain  the  absolute  best 
parameters  is  not  attempted  because  there  will  be  some 
variation  in  the  data  from  mission  to  mission. 

The  next  step  in  developing  the  model  is  the  calculation  of 
the  single-T  power  spectrum.  Although  we  have  chosen  specific 
parameters  to  represent  the  final  spectrum,  we  will  still 
investigate  the  shapes  of  F«p(f)  for  the  same  q,  s  and  p  as  in 
figures  13a-c  and  l^a-c.  Figures  l6a-c  are  the  corresponding 
normalized  spectra  FT(f)  versus  normalized  frequency  f/f^ 

The  effect  of  power  in  the  spatial  spectrum  below  K-=2iTf/v  is 
ignored  because  of  the  lower  limit  of  the  integral  ( 1^) .  Thus 
there  is  net  as  pronounced  a  difference  between  the  temporal 
spectra  as  there  was  w ith  the  spatial.  There  is  sufficient 
low  frequency  power  however  to  follow  the  single-T  spectral 
data  of  mission  210873-  Figure  17  shows  that  comparison,  and 
it  is  apparent  how  the  parametric  values  q=2  and  s-1  have  caused 
a  more  gradual  break  into  the  low  frequency  asymptote. 

To  round  out  the  comparison  between  the  new  model  and  the 
data  of  mission  210873,  we  present  the  temperature-difference 
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Figure  17.  Single  probe  temperature  spectral  data  for 
mission  210373  and  theoretical  curve  based  on  the  model 
for  q=2,  p=11/6,  and  s=1. 


spectra  in  figure  18.  The  measured  average  orientation  angle 
6  of  55.6°  is  used.  There  is  very  little  noticeable  difference 
between  these  curves  and  those  in  figure  which  were  for  the 
von  Karman  spectrum.  This  is  as  expected  since  spatial 
spectrum  changes  below  k=L0_1  are  not  noticed  by  a  probe 
separation  r<<LQ.  The  remaining  error  must  be  due  to  the  use 
of  Taylor's  hypothesis. 

Shkarovsky's  Model 

In  a  paper  on  the  generalized  spatial  spectrum, 

Shkarovsky [22]  develops  a  spectrum  which  has  two  power  law 
dependencies,  such  as  used  here.  It  was  not  the  basis  for 
the  current  analyses  since  it  lacks  a  parameter  to  adjust  the 
shape  of  the  transition  region,  fitter  developing  the  proper¬ 
ties  of  this  familiar  spectrum: 

VA<)  =  [1+(K/Ki)2]-U/2.,  (57) 

he  multiplies  two  such  spectra  which  have  k  <<ic  : 

V(k)  =  [1+(k/k  )2]“v/2[1+(ic/k  )2](m_ij  )/2  ,  (58) 

1  2 

In  the  very  lowest  wavenumbers  this  model  approaches  unity, 
a  property  which  is  not  necessary.  We  eliminate  that  low 
wavenumber  asymptote  by  letting  k»k  : 

VU)=c  k_u[1+(k/<  )2]^‘-w')/2.  (59) 

1  2 
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The  two  parameters  are  p  and  p^,  the  power  law  dependencies  in 

the  low  and  high  frequency  regions,  respectively.  We  can 

introduce  a  third  parameter  t  to  adjust  the  transition  region: 

V(ic)=c  <"p[1+(kA  )t](p_p")/t.  (60) 

1  2 

This  can  now  be  compared  with  the  model  of  this  paper,  equation 
(i<5),  rewritten  as 

*t(k)=  C2K-sP[l+(KL0)q"S]-P.  (61) 

The  comparison  of  (60)  with  (61)  is  obvious: 

<  *L  ,  c  =c 
2  o  1  2 

t=q-s 

p^qp 

p=*sp  and 

p*(p'-p)/t 

s=pt/(p'-p) 

q*P '  t/(p'-p) .  (62) 

If  we  take  the  optimum  parameters  q=2,  s=l  and  p=ll/6  as  being 
correct,  then  p=ll/6,  p'=ll/3  and  t=l.  Unfortunately, 
Shkarovsky’s  model  has  t  fixed  at  2  so  his  model  is  not 
reducible  to  the  optimum  form.  Alternatively  if  we  assume  his 
t=2  is  correct  and  furthermore  use  p'*ll/3  and  ps'llA  then 
p=ll/12,  q«4  and  s*2.  The  q*1!  tells  us  that  the  knee  of  his 
curve  is  considerably  more  abrupt  than  found  acceptable. 


VI.  COMPARISON  OF  MODELS  WITH  OTHER  DATA 


Since  one  mission  at  one  site  is  not  sufficient  to  develop 
a  generalized  model,  other  data  will  now  be  compared  with  the 
new  theory  to  show  an  equally  good  fit.  The  missions  are: 

200373  at  Kirtland  AFB  where  sensors  were  at  20m  and  33m  above 
uniform  terrain,  and  0^0973  and  011073  at  Griffiss  AFB  under 
ideal  conditions  similar  to  210873* 

Mission  200373 

To  see  if  local  terrain  governed  the  low  frequency  behavior 
of  the  spectra  at  Griffiss,  data  taken  at  Kirtland  AFB  were 
processed  for  spectral  information  only.  Two  sensor  pairs 
were  positioned  at  20m  and  33m  above  ground  on  a  stable  tower. 
The  tower  is  near  the  center  of  a  shallow  canyon  1.6km  wide 
and  with  walls  ^.m  and  27*m  above  the  base  of  the  tower. 

There  is  sparse  desert  vegetation  -  low  brush  and  juniper  trees. 
The  canyon  floor  is  hard  dirt  and  is  rocky.  Winds  are 
typically  up  or  down  canyon  and  the  ratio  ov/v  is  typically 
low.  Whereas  the  nine  sensors  at  Griffiss  are  on  a  long 
manually  rotatable  boom,  the  sensor  pairs  at  Kirtland  are  on 
windvanes  which  track  direction  variations  at  about  0.5  Hz  with 
a  threshold  speed  of  0.5  m/sec.  The  mission  summary  is 
in  figure  19  for  only  the  33m  altitude. 


To  save  space  in  the  spectral  plots,  ecch  page  will 
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Figure  19. 

Mission  Summary  200373 
(Values  for  33m  altitude) 

Cloud  Cover 

45KFT  Scattered 
25KFT  Scattered 

Visibility  30mi. 

Barometric  Pressure  811  mb 

Dev  Point  4.4°C 

Temperature  15.6°C 

Ground  Conditions*  Dry 
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represent  a  single  spectrum  for  both  altitudes.  However,  the 
theoretical  curves  will  apply  only  to  the  33m  altitude  data. 

In  short  these  curves,  figures  20a-f ,  offer  no  surprises  but 
corroborate  the  use  of  the  new  model.  The  temperature-differ¬ 
ence  spectrum  in  figure  20c  still  has  the  up-turn  in  the  data 
below  the  external  scale  frequency  f ^ - 

Mission  040973 

Unfortunately  the  Kirtland  mission  could  rot  include  a 
measurement  of  the  structure  function,  so  tj  get  further  back¬ 
up  spatially,  another  mission  was  run  at  Griffiss  AFB  for 
structure  function  data  only.  The  mission  summary  in  figure 
21  shows  that  occasionally  the  wind  came  along  the  line  of 
sensors.  Those  sections  of  data  were  not  averaged  into  the 
overall  structure  function  and  in  the  summary  are  marked  with 
astericks.  The  normalized  structure  function,  measured  with 
only  seven  sensors  operative,  is  plotted  in  figure  22.  The 
fit  of  the  data  to  the  modeled  curve  is  certainly  commensurate 
with  the  fit  of  mission  210873- 

> 

Mission  011073 

The  most  nearly  ideal  mission  is  saved  for  last.  This 
mission  at  Griffiss  had  low  wind  soeed  and  direction  variances 
and  thus  we  expect  the  frozen-flow  hypothesis  to  work  well. 

The  summary,  figure  23,  shows  this.  The  structure  function 
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Figure  20a*  Single  probe  spectral  data  for  mission  200373 
Kirtland  AF3,  20m  altitude,  and  a  theoretical  curve  based 
on  the  proposed  model. 


Figure  20d.  Modulus  of  the  temperature  cross-spectrum,  data 
and  theory  based  on  the  proposed  model. 


log(f) 

figure  20e.  Phase  component  of  the  cross- spectral  data.  The 
deviation  of  these  data  from  zero  for  f?4  Hz  demonstrates  a 
slight  misalignment  of  the  probes  on  the  wind  vane. 


Figure  21 . 

Mission  Summary  040973 
Cloud  Cover  0# 

Visibility  30mi. 

Barometric  Pressure  1005  mb 
Dew  Point  19.7°C 
Temperature  32°C 
Ground  Conditions:  Dry 
Sensor  Height  2.26m 
Grass  Height  0.07m 

v*»1 .6  2m/ sec 
ov=0. 488m/ sec 
5*1 15.6  °s»2 . 02rad 
og=*0c43rad 
AQ=2.92m 
C==0.142°C3m"a/3 
o^/ya 0.30 

r=0,243m 


♦Because  of  high  wind 
direction  variance, 
spectra  for  these  times 
were  not  used. 
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T 
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Figure  22.  Normalized  temperature  structure  function  data  and 
theoretical  curve  based  on  the  proposed  model.  Mission  040973 
at  Griffiss  AFB. 
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Figure  23. 

Mission  Summary  011073 

Cloud  Cover  20$, 

Altitude  4000ft. 

Visibility  30mi. 

Barometric  Pressure  1028  mb 

Dew  Point  7.2°C 

Temperature  22cC 

Ground  Conditions:  Dry 

Sensor  Height  2.26m 

Grass  Height  0.07m 

v=*3.16m/sec 
cr^aO. 775m/  sec 
B»950«1.66rad 
a0»O.33rad 
Ao»5.11m 
C|»0.26°C3m"a'3 
ct^Tv-0.25 


r*0.428m 


and  spectral  data  in  figures  2*Ja-f  need  little  elaboration 
as  they  show  the  same  behavior  as  previously  presented  missions. 
One  new  aspect  is  that  the  difference  spectrum  below  the 
external-scale  frequency  shows  an  excellent  fit  to  the 
theoretical  curve.  This  is  a  result  of  low  values  of  oy/v 
and  Og. 


y 
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Figure  24a.  I^rmalized  temperature  structure  function  data 
and  theoretical  vvrve  based  on  the  nroposed  model.  Mission 
011073,  Griffis s  APB. 
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Figure  24b.  Single  probe  temperature  spectral  data  and  theore 
tioal  curve  based  on  the  proposed  model.  Mission  011075  at 
Griff iss  APB. 
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log(f) 

Figure  24c.  Temperature  difference  spectral  data  and  theoreti¬ 
cal  curve  based  on  the  proposed  model. 


VII.  SUMMARY  OP  EQUATIONS  FOR  THE  NEW  MODEL 


This  section  summarizes  the  equations  for  the  new 
parameters  q=2,  s=l,  p=ll/6.  The  new  spatial  spectrum  is 

*TU)  =  0.0330  C2  (<2+k/L0)“1i/6.  (63) 

The  structure  function  is  calculated  by  substitution  of  (63) 
into  (8): 

DT(r)=cf,  rl/3[l.l078(£~) 

Im  u(f  ,  3  ,  l^)] 

DT(r<<Lo)=0|  r2//a, 

DT(Lo«r)=1.1078  Lq2/jC|, 
and  DT(r«Lo)/DT(L0«r)  =  0.9027(^-) 

The  function  U(a,b,c)  is  a  confluent  hypergeometric  function[15] . 
The  external  scale  A0  in  terms  of  L0  is 

A0  =  1.1660  Lq  (65) 

and  is  not  much  different  from  A0  *  1.071Lo  for  the  von 
Karman  spectrum. 

The  single-T  power  spectrum  does  not  have  as  neat  a  form 
as  the  structure  function  and  is  best  obtained  by  numerical 
integration  of  (63)  in  (14).  Formally  it  may  be  represented  as 

PT(f)-0. 07307  C*.  v2/j  fs/3  F  (11/6,  5/3;  8/3;  -f  /f) 

1  21  1 

(66) 


_ii-6i73(y 

(644) 


2/3 
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where  F  (a,b;c;z)  is  the  hypergeometric  function[15] .  The 
21 

two  asymptotes  for  FT(f)  are 

liir,  Ft( f )=13 .919  v-1L05/3[l-1.1232(f  /f)~l/s]  (67) 

f-K)  1 

and 

iim  FT(f)=0. 07307  C,|  v2/s  f-5'''3.  (68) 

f-*» 

The  inertial  subrange  dependence  (68)  is  precisely  the  same 
as  the  dependence  for  the  von  Karman  spectrum  (15),  as  it  must 
be . 


The  temperature-difference  spectrum  can  only  be  obtained 

by  numerical  integration  of  (63)  in  (20).  The  two  asymptotes 

f<<f  ,  f  and  f>>f  ,  f  are  calculable.  When  f>>f  the 
1212  1 

modified  portion  of  the  spatial  spectrum  has  been  ignored  and 

the  form  derived  from  the  von  Karman  spectrum  is  applicable. 

This  is  equation  (2*1)  except  now  in  (26)  3  is  Riven  by 

(fsin9)/f  .  The  asymptote  for  small  f  where  also  6^0  and 
2 

rsin 0  <<L0  is  found  via  Erdelyi  et  al[23],  transform  pair  10  on 
page  23: 

Iim  FAT(  f )-5 . 8276  v-1C^.  (rsinO)5/3 

f’°  [1-1. M5  ].  (69) 

This  assures  us  that  the  new  modifications  to  the  theory  do 
not  account  for  a  excess  low  frequency  power  in  the  difference 
spectrum,  since  this  is  identical  to  (2*1)  for  the  same 


conditions . 


VIII.  CONCLUSIONS 


In  this  paper  we  have  developed  a n  empirical  model  for  the 
spatial  temperature  spectrum.  The  spatial  spectrum  was  not 
measured  directly,  rather  the  spatial  structure  function,  which 
is  related  via  a  Fourier  transform  to  the  spectrum,  was 
measured.  A  spectral  model  was  proposed  and  the  parameters 
varied  until  a  best  fit  to  the  structure  function  data  was 
obtained.  The  best  model  was  then  transformed  to  the  temporal 
spectrum  using  the  frozen  flow  hypothesis  and  again  good 
agreement  was  found.  The  only  residual  error  traceable  to 
inappropriateness  of  frozen  flow  was  in  the  difference  spectrum 
and  e.ven  that  error  was  reduced  when  the  variances  of  wind 
direction  and  speed  lessened. 

Four  data-gathering  missions  were  run  at  the  Kirtland  AFB 
and  Griffiss  AFB  laser  propagat  sites,  both  on  rather 
uniform  terrain.  Even  for  substantially  different  conditions 
at  these  two  sites,  tne  data  fit  the  new  model  well.  An 
Inspection  of  the  structure  function  data  at  Griffiss  revealed 
no  not j cable  dependence  of  curve  shape  or  external  scale  on 
wind  direction  wi’.b  respect  to  the  sensor  line.  Presently 
there  is  no  proper  experiment  to  measure  the  three-dimensional 
external  scale  since  sensor  wake  effects  are  significant 
when  a  measurement  along  the  wind  vector  is  attempted  and  since 
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the  external  scale  varies  with  altitude  in  a  vertical 
measurement. 

All  the  data  exhibited  the  Kolmogorov  behavior  in  the 
inertial  ranges.  This  is  a  result  of  running  under  the  best 
conditions  for  the  available  sites:  moderate  winds,  well- 
developed  turbulence  and  convection  at  mid-day,  good  mixing, 
few  clouds  and  low  wind  direction  fluctuations.  A  different 
behavior  is  expected  at  the  near-neutral  conditions  near  dawn 
and  dusk  and  in  the  stable  conditions  usually  encountered  at 
night.  Measurements  made  by  Kaimal  et  al[24]  in  a  Kansas 
wheat  field  indicate  the  more  nearly  stable  conditions  may  fit 
the  von  Karman  model  rather  than  the  proposed  one.  However, 
in  unstable  conditions  the  data  seem  to  support  the  new  model. 
Further  experimentation  is  needed  to  determine  if  the  low- 
frequency  behavior  depends  on  the  Richardson's  number  or  on 
the  temperature  structure  constant.  Ur:der  some  conditions 
the  turbulence  may  be  non-Kolmogoro?.  It  is  hoped  the 
generality  in  the  early  development  of  the  model  will  allow  the 
reader  to  extend  this  model  to  his  particular  case  if  necessary. 

The  new  model  seems  to  make  the  mathematics  more 
difficult,  but  it  does  not  introduce  new  mathematical 
difficulties  such  as  non-convergence.  Those  propagation 
calculations  which  depend  on  large  scales  of  turbulence,  such 


as  beam  wander,  should  be  redone  using  the  new  model.  If  the 
calculations  do  not  depend  on  large  scales  but  require  a 
formula  which  has  finite  power,  then  the  von  Karman  spectrum 
should  be  used  as  it  is  easier  to  work  with. 
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